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ABSTRACT 



^ Aims. We study the CO and the radiocontinuum emission in an active galaxy to analyze the interplay between the central activity and the 
molecular gas. 

Methods. We present new high-resolution observations of the '"CO(l-O) and '"C0(2-l) emission lines, and 3.5 cm and 20 cm radio continuum 
] emission in the central region of the LINER/starburst galaxy NGC 6764. 

. • Results. The galaxy has an outflow morphology in radio continuum, spatially coincident with the CO and Ho- emission, and centered slightly 
C ' off the radio continuum peak at the LINER nucleus. The total molecular gas mass in the center is about 7x10^ Me, using a CO luminosity to 
total molecular gas conversion factor that is three times lower than the standard one. '^CO(l-O) emission is found near the boundaries of the 
radio continuum emission cone. The outflow has a projected expansion velocity of 25 km s 'relative to the systemic velocity of NGC 6764. 
■ About 4 X 10* Mo of molecular gas is detected in the outflow. The approximate location (~ 1 kpc) of the dynamical inner Lindblad resonance 
^ , has been derived from the rotation curve. The peak of the CO emission is slightly (< 200 pc) offset from the peak of the radio continuum. 
! '~j Conclusions. The molecular gas has most likely been ejected by the stellar winds from the recent starburst, but the CO line ratios show 
K> indication of an interaction with the AGN. The energy released by the nuclear starburst is sufficient to explain the observed outflow, even if the 
^ data cannot exclude the AGN from being the major energy source. Comparison of the outflow with hydrodynamical simulations suggests that 
■ - - ' the nuclear starburst is 3-7 Myr old and the bubble-like outflow is still confined and not freely expanding. 

Key words. ISM: jets and outflows — Galaxies: active — Galaxies: ISM — Galaxies Individual: NGC 6764 — Galaxies: starburst — Radio 
continuum: galaxies 



1. Introduction 

The nuclear region surrounding active galactic nuclei (AGN) 
often has the highest level of molecular gas concentration and 
star formation in galaxies. The processes that drive the excita- 
tion and kinematics of the circumnuclear gas, namely, the abun- 
dant star formation and the jets and photoionisation associated 
with the AGN, influence the circumnuclear gas in several ways. 
The onset and shutdown of the starburst are driven by the prop- 
erties of the molecular gas, including its density and velocity 
dispersion, by gravitational instabilities (e.g.. Combes 2001), 
and by processes that tend to prevent star formation, such as 
the local tidal field (Launhardt, Zylka & Mezger 2002), star- 
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burst feedback (e.g., Dopita 1985, Scalo & Chapell 1999), gas 
consumption (e.g., Franceschini et al. 1998), and tidal shear 
(e.g., Kenney et al. 1993). Fueling of starbursts is facilitated 
by nonaxisymmetric components such as bars, and by tidal 
torques from interactions and mergers (e.g. Combes & Gerin 
1985, Athanassoula 1992, Mihos & Hernquist 1994, Laine et 
al. 2002) 

The formation of powerful gas and plasma outflows from 
the center of disk galaxies provides the basis for the metal 
enrichment of galaxy halos (e.g., Kunth et al. 2002). Several 
mechanisms have been proposed to explain these outflows. 
The first possibility is a thermal wind from a circumnuclear 
starburst (e.g., Colina, Arribas & Borne 1999). The second 
possibility is a thermal wind from an AGN, as suggested by, 
e.g., Krolik & Bergelman (1986). Finally, it is possible that 
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ram-pressure along the radio jet is driving the outflow (e.g., 
Taylor, Dyson & Axon 1992). Large-scale radio halos out- 
side the galaxy disk, which most likely result from plasma and 
magnetic fields that have been blown off" the disk by one of 
the above-mentioned mechanisms, have been known for more 
than two decades (e.g.. Hummel, van Gorkom & Kotanyi 1983, 
Hummel, Beck & Dettmar 1991). 

The nuclear starburst leads to the formation of OB asso- 
ciations or super star clusters, from which the massive stars, 
through strong supersonic stellar winds, inject both energy and 
mass into the surrounding interstellar medium. The formation 
of an expanding starburst "superbubble" will be accompanied 
by extended X-ray emission from hot gas, and optical Ho- line 
emission from the cooler shell. The superbubble that has a 
larger internal ISM pressure than its surroundings expands into 
the low-pressure halo ISM until the onset of Rayleigh-Taylor 
instabiUties, which disrupt the superbubble and allow the hot 
gas to escape in the form of a galactic wind, as is seen, e.g., in 
M82 (Strickland, Ponman & Stevens 1997). The cold molec- 
ular gas is expected to be compressed around the superbubble 
into an expanding cocoon. 

Molecular gas in outflows has been detected in a few galax- 
ies. M82 has an outflow bubble expanding at 45 km s"'(Wei6 
et al. 1999), with a molecular outflow gas mass of ~ lO*" Mq. 
NGC 2782 harbors an outflow with about ~ 2 x 10^ AIq of 
molecular gas (Jogee et al. 1998). The outflow in NGC 2782 
appears to be in the early stages (its age is ~ 4 x 10^ years), so 
that it is seen before the onset of Rayleigh-Taylor instabilities. 
In NGC 3079 four expanding CO shells have been identified 
(Irwin & Sofue 1996). These appear to have originated from a 
single starburst that took place a few Myrs ago. 

The behavior of the molecular gas in such energetic out- 
flows is still poorly understood. The kinematic and potential 
energies of the molecular gas are part of the total energy bud- 
get of the outflow, and one can use the observed properties of 
the gas (its location, mass and velocity) to constrain the ori- 
gin of the outflow. Recently Sakamoto et al. (2006) observed 
molecular gas shells with the Submillimeter Array (SMA) in 
the circumnuclear disk of NGC 253. Further constraints on the 
nature and location of the energy source can be provided by 
the current gas phase, as in NGC 2782 (Jogee et al. 1998), or 
an unusual geometry of the outflow, such as the off-center out- 
flow in NGC 5775 (Lee et al. 2001). Since gas is fuel for star 
formation, when it falls back onto the galaxy disk it may induce 
new episodes of star formation (e.g., Norman & Ikeuchi 1989). 

NGC 6764 is a barred spiral with an optically classi- 
fied low-ionisation nuclear emission-line region (LINER) nu- 
cleus (Osterbrock & Cohen 1982). This nearby active galaxy 
(D=32 Mpc; 1" - 160 pc) hosts a very young nuclear star- 
burst as shown by its Wolf-Rayet (W-R) feature at 0.47 fim 
(Hell; Schinnerer et al. 2000). Eckart et al. (1991, 1996) and 
Schinnerer et al. (2000) analyzed the central starburst, traced 
by the W-R emission, in great detail. The millimetre obser- 
vations have revealed a central concentration of molecular gas 
which also emits near-infrared line emission. The 2.12 yum H2 
line has features that are also seen in radio continuum. A pop- 
ulation synthesis model based on NIR spectral lines suggests 
that NGC 6764 has undergone two recent starbursts. One of 



Table 1. Global properties of NGC 6764. 



Property 


Value 


References 


Hubble type 


SB(s)bc 


NED 


Distance (Mpc) 


32 


Ho = 75 km s"' Mpc"' 


Inclination 


61.5° 


LEDA 


Major Axis P.A. 


62° 


LEDA 


Bar P.A. 


74.6° 


LEDA 


Si2;um (Jy) 


0.36 


NED 


S25yum (Jy) 


1.29 


NED 


SfiOyum (Jy) 


6.33 


NED 


Siooyum (Jy) 


11.56 


NED 


LpiRdO'" Lo) 


1.1 




Lh„ (10* erg s-') 


23 




Mhi (lO'Mo) 


7.0 


NED 


S900,, (mJy) 


309 


WENSS 


Ssocm (mJy) 


244 


Douglas et al. (1996) 


S20cm (mJy) 


133 


White & Becker (1992) 


Seem (mJy) 


34 


Becker, White & Edwards (1991) 


S3.5<.m (mJy) 


31 


This paper 



them took place 3-5 Myr ago, producing W-R features, and 
another 15-50 Myr ago. RO SAT X-iay data trace mainly the 
central AGN which varies by a factor of two over seven days 
(Schinnerer et al. 2000). By studying HCN and CO emission 
in the center of NGC 6764, Kohno et al. (2001) classified 
this galaxy as a combined AGN/starburst, based on the high 
HCN/CO line ratio (~ 0.7) in the nucleus which suggests the 
presence of nuclear activity. Recently Hota & Saikia (2006, 
hereafter HS06) presented GMRT+VLA radiocontinuum and 
HI observations together with Ha data of the galaxy NGC 6764 
where they found bipolar radio bubbles along the minor axis to- 
gether with a compact core. Their HI observations reveal two 
peaks of emission at the ends of the bar and a depletion in the 
center together an HI absorption towards the radio core. The 
coexistence of an AGN and a starburst nucleus in NGC 6764 
makes this galaxy an ideal target for studying the influence of 
both processes. 

We present new high spatial resolution observations of the 
'2CO(1-0) and '2cO(2-l) lines made with the IRAM Plateau 
de Bure Interferometer (PdBfl). We also present high resolu- 
tion vlaI radio continuum maps at 3.5 cm and 20 cm (Laine 
et al. 2006). We describe our observations in §2. In §3 we dis- 
cuss the thermal and synchrotron emission associated with the 
hot ISM, and their connection to the nuclear star formation. 
The distribution, kinematics, physical properties, and mass of 
the molecular gas in the outflow in NGC 6764 are described in 
§4. The energies involved in the outflow as determined by the 
molecular gas and radio continuum observations of the outflow 
are addressed in §5. We discuss our results in §6 and summarize 
them in §7. Throughout this paper, we use the following cos- 
mological parameters: Hq = 75 km s"', = 0.7 and Q;, = 0.3 

' Based on observations carried out with the IRAM Plateau de 
Bure Interferometer. IRAM is operated by CNRS (France), the MPG 
(Germany) and the IGN (Spain). 

- The VLA is a research facility of the National Radio Astronomy 
Observatory (NRAO), which is operated by Associated Universities 
Inc., under contract with the National Science Foundation. 
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2. Observations and data reduction 

2.1. Radio continuum data 

The radio continuum data were obtained with the Very Large 
Array (VLA) at 3.5 and 20 cm in C and A configurations 
(see Figure [T]), respectively. We obtained a comparable spa- 
tial resolution at the two observed wavelengths, r.'39 x 1'.'26 
(P.A.=-45°) at 20 cm and 2'.'50 x 2" (P.A.=15.5°) at 3.5 cm. 
For further details about the radio continuum observations, see 
Laine et al. (2006). The peak position of the 20 cm emission, 
a20cm(^2000) = 19:08:16.316, 520cm(^2000) = 50:55:59.37, 
was determined to an accuracy better than 0'.' 1 , by filtering out 
large spatial frequencies, using a wavelet transformation (e.g., 
Leon et al. 2000), and fitting a Gaussian to the remaining emis- 
sion. 



2.2. CO data 

The '2co(l-0)and '2co(2-l) lines were observed with 
the IRAM interferometer (see Table |2|i, located on the 
Plateau de Bure (France), in the B and C configurations 
of the array, in 1998. We calibrated the visibilities us- 
ing the radio galaxy 3C 380 as an amplitude calibra- 
tor and 1823-1-568 and 3C 380 as phase references. The 
Plateau de Bure Interferometer (PdBI) was phase-centered at 
ay2000=19:08: 16.370 and 5y2000=+50:55:59.58. The velocity 
resolution was set at 10 km s ' for the two CO transitions, 
and the CLEANed maps have been restored with a synthesized 
beam of 2'.'0 x 1'.'6 [PA.=7°; '200(1-0)], and 1'.'25 x 1'.'25 
[P.A.=0°; '^€0(2-1) ]. The rms noise for the channel maps is 
2 mJy beam-i [58 mK; '^00(1-0) ], and 5.5 mJy beam"' [82.5 
mK; '^€0(2-1) ]. Hereafter the velocity reference is fixed to 
2400 km s"', giving a redshifted frequency of 114.3489 GHz 
for the i2cO( 1-0) observations, and 228.6932 GHz for the 
'2CO(2-l) observations. 

The total integrated '^€0(1-0) flux is 225 Jy km s^Hn the 
55" beam of the NRAO 12-m telescope (Sanders & Mirabel 
1985). Our PdBI '^00(1-0) observations detected a total flux 
of 112 Jy kms"', i.e., filtering out about 50% of the total 
flux, assuming all the emission is concentrated in the 45" pri- 
mary beam of the PdBI. This implies that a large fraction 
of the molecular gas must be in an extended smooth compo- 
nent to which the interferometer is insensitive. From Eckart 
et al. (1991) we estimate with IRAM-30m observations the 
^^CO{2-\) integrated intensity to be ~ 480 Jy km s 'in the cen- 
ter of NGC 6764, which implies a filtering of about 65% of the 
'^€0(2-1) flux by the PdBI observations. 

The integrated intensity maps for the '^€0(1-0) and 
'2CO(2-l) emission lines are shown in Figure |2] The posi- 
tion of the CO emission peak q'co(^2000) = 19:08:16.402, 
5co(-/2000) - 50:56:0.29 has been estimated using the higher 
spatial resolution of the '^00(2-1) transition. Assuming that 
the spiral arms are trailing, and noting that the receding ve- 
locities in the disk are on the eastern side of NGC 6764, the 
northern side of the galaxy disk is closer to us. 



Table 2. PdBI observations 



Parameter 


'=CO(1-0) 


'2CO(2-l) 


Field of View 


45"(6.9kpc) 


22"(3.4 kpc) 


Bandwidth! 


270 MHz 


473 MHz 


Velocity resolution 


lOkms' 


lOkms"' 


Syntliesized beam size 


2" x 1'.'6, PA. T 


1'.'25 x 1'.'25 


rms noise (per channel) 


2 mJy beam"' 


6 mJy beam"' 



3. Thermal and synchrotron emission 

3.1. Radio continuum 

The total integrated 20 cm flux density from the VLA obser- 
vations is 90 mJy, which compared to the single dish flux den- 
sity of 133 mJy (White & Becker 1992), impHes that -70% 
of the total flux density is recovered by the VLA observations. 
The total 3.5 cm flux density is 31 mJy. Both VLA maps at 
20 and 3.5 cm present the same features with a better spatial 
resolution at 20 cm (see Fig.[T]). A bright peak is observed to- 
wards the nucleus of NGC 6764. Two 5" spurs emerge from 
the unresolved nucleus, towards the south and the south-east. 
The emission near the nucleus is not consistent with a just a 
nuclear point source but hints at the presence of a nuclear radio 
jet, as proposed by HS06. The most remarkable feature is the 
shell-like emission, at both wavelengths, centered on the cen- 
tral peak with its major axis perpendicular to the major axis of 
NGC 6764 and extending about 1 kpc above and below the disk 
plane of the highly inclined galaxy NGC 6764. The clumpy 
shell emission is brighter towards south than towards north. 
These radio bubbles were previously reported by Baum et al. 
(1993), HS06 and Laine et al. (2006). 

3.2. Ha emission 

An Ha image of NGC 6764 was kindly given to us by D. 
Frayer. The observations were done at the Palomar 60-inch 
telescope with a seeing of 1-3", using the X-filter centered at 
6616 Awith a width of 20 A, leaving out most of the [Nil] 
emission line. Figure [3] shows the Ho- emission at the center 
of NGC 6764. Apart from the most intense Ha emission in 
the central region, two other regions of emission are seen near 
the tip of the bar, which also lies near the beginning of the 
spiral arms. The most striking characteristic of the Ha emis- 
sion near the center of NGC 6764 is its bubble-like shape, but 
its size is smaller than the size of the radio continuum out- 
flow feature in the south. In Figure |4] a 4" spur is observed 
to the southeast of the location of the '2CO(2-l) peak, bend- 
ing towards south. This spur is spatially correlated with faint 
(< 6 Jy beam"' km s"') '2CO(2-l) emission. 

In Figure |5] the Ho- emission is compared with the radio 
continuum spectral index. The "hot spot" in the bubble cavity 
with a flat spectral index (a ~ -0.2) is associated with Ha 
emission driven by the nuclear starburst. Similarly, the tip of 
the southeastern Ha emission spur coincides with a small "hot 
spot" that has a flat spectrum. 
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19''0B"17.0= 16.5= 16.0' 15.5= 19''0B='17.0= 16.5= 16.0' 15,5= 

Right Ascension (J2000) Riglit Ascension (J3000) 

Fig. 1. Left: radio continuum emission at 20 cm with a synthesized beam of 1'.'39 X 1'.'26. The contour levels are at (0.1, 0.5, 
1, 1.5, 2, 3, 4, 8, 12) mJy beam"' and the rms noise is 70 yuJy beam""'. Right: radio continuum emission at 3.5 cm (VLA) with 
a synthesized beam of 2'.'50 x 2'.'0. The contour levels are at (0.1, 0.2, 0.4, 0.8, 1, 2, 4, 8) mJy beam"' and the rms noise is 56 
fiJy beam"'. 




19'0B"17- IB- 15- ig'OB-lV.O' 16.5' 16.0= 

Right Ascension (IZOOQ) Kight Ascension (J2000) 



Fig. 2. Left: VLA 20 cm radio continuum emission (greyscale) overlaid on the integrated intensity image of the 

'^CO(l-O) emission. Right: '^€0(2-1) integrated intensity image. The contours for '2CO(1-0) are at (1, 2, 4, 8, 15, 20, 25) 
Jy beam ' km s ', and for '^CO(2-l) at (2, 4, 6, 12, 24, 52,) Jy beam ' km s '. The position of the major axis of the galaxy and 
of the bar are shown by soUd and dashed Unes, respectively. 



3.3. Nuclear star formation 

The nuclear star formation history was studied in detail by 
Schinnerer et al. (2000) using near-infrared K-band imaging 
spectroscopy of the very center of NGC 6764. To get the cur- 



rent star formation rate we have used here the FIR and the Ha 
emission line. We estimate as well the contribution of the ther- 
mal and non-thermal radiocontiuum emission. 
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+ 50°57'00" - 




19''08""20^ IS'' 
Right Ascension (J2000) 



Fig. 3. B image (from Digital Sky Survey) of NGC 6764 over- 
laid by Ha emission contours. The contours are spaced by 0.5 
dex steps, starting at the 0.5 cr level. 



+50°56'10" - 




°17.5' 17.0= 16.5' 16.0' 15.5' 15.0' 

Right Ascension (J2000) 

Fig. 4. Ha emission (greyscale) overlaid by the 
'^CO(2-l) emission contours (same levels as in Fig.|2]i. 



3.3.1. Star formation rate 

The total FIR luminosity is computed for the whole galaxy 
from the 60 and 100 yum IRAS fluxes (see Table [T]), but is very 
close to the bolometric value Lboi - 9 x 10^ L© computed 
for the nuclear region by Schinnerer et al (2000). Following 
Condon (1992), the star formation rate (SFR) for massive stars 
is. SFR{M > 5 Mq) = 9.1 X \Q-^\Lf,rILq), which gives a 
SFR of 1 Mo/yr 




19''0B""17.0' 16.5' 16.0' 

Right Ascension (J2000) 

Fig. 5. Spectral index -a between 3.5 cm and 20 cm 
(greyscale), overlaid by Ha emission contours. 

Table 3. Outflow bubble and nuclear contributions to the Ha 
and 20 cm radio continuum emission derived by fitting a 2D 
Gaussian to the central emission. The molecular gas mass for 
the NSB is the mass encompassed by one '^CO(l-O) beam, and 
has a projected size of 3 10x250 pc. 





Ha 


20 cm 3.5 cm 


M(H2) 




(10'' erg.s-') 


(mjy) (mjy) 


(VHo) 


Nucleus 


~ 130 


90 28 


4.3 x lO** 


Outflow 


~ 45 


45 12 


1.7 X 10' 



The Ha line is mainly emitted by HII regions where newborn 
massive stars are ionizing the gas. The total Ha flux can be 
used as a tracer of the star formation rate (Kennicutt 1998): 



L(Ha)(erg/i) = 1.26 x Kf^SFRiMelyr) 



(1) 



The Ha flux is 1.3 10"*' erg/s for the nuclear starburst for an 
aperture size of 3"(see Table|3]l, which gives a SFR of 1 Alo/yr, 
consistent with the SFR derived from the FIR luminosity. These 
two derivations of the SFR take into account only massive 
stars. Including all stellar masses, and accounting for the ef- 
fects of optical extinction, Eckart et al. (1991) derived a SFR 
of 4 Alo/yr for the nuclear starburst (NSB), and a supernova 
rate of 1.2 x lO^^/yr. 

3.3.2. (Non-)Thermal radiocontinuum 

While the radio continuum emission at long wavelengths is 
mainly produced by the compact source (AGN) the contribu- 
tion from the NSB is important at 3.5 cm creating a flatter 
spectrum. By fitting the radio continuum between 6 cm and 80 
cm where the non-thermal contribution is dominant, we found 
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Table 4. Thermal (free-free) and non-thermal (synchrotron) 
contribution to the observed radio continuum emission at dif- 
ferent wavelengths at the center of NGC 6764. The thermal 
contribution is estimated from the Ha luminosity (see text). 



Wavelength 


Thermal 


Non-thermal 


ratio 


Thermal 
Total 




(mJy) 


(mJy) 




20 cm 


17.6 


115.4 


13% 




6 cm 


15.6 


18.4 


46% 




3.5 cm 


14.8 


16.3 


48% 





1 ,00 




> 0.70 r 

j| ; 

l/f - 
0.60 r 



0.50 f - 

5 10 15 20 25 30 

A (cm) 

Fig. 6. Ratio of the thermal radio continuum emission (mJy) 
and Ha luminosity (10"^*^ erg/s) at the distance of NGC 6764, 
assuming an electron temperature of 10,000 K and an optical 
extinction of Ay -4. 



a spectral index of a = -0.46. Extrapolating that fit down to 
3.5 cm, we obtained a non-thermal flux density of 6 mJy, which 
is much below the measured value of 3 1 mJy. 

The free-free contribution ("bremsstrahlung") can be com- 
puted from the Ha flux. The Balmer lines provide an estimate 
of the number of ionizing photons present in the source. If the 
optical extinction is known, the Ha flux can be used to compute 
the thermal radio continuum flux (e.g., Caplan & Deharveng 
1986). Assuming an electron temperature of lO'* K and an opti- 
cal extinction of Ay - 4 (Eckart et al. 1991), the thermal radio 
continuum flux density at 3.5 cm is 14.8 mJy, obtained by in- 
tegrating the Ha luminosity in the outflow. The Ho' extinction 
is given by Afja - 0.605 * Ay- Figure |6] shows the ratio be- 
tween the radio flux density (mJy) and the Ha luminosity (in 
units of 10"^*^ erg/s) for the thermal (free-free) contribution ver- 
sus wavelength. This ratio is 75% at 3.5 cm. The maximum 
thermal contribution to the flux density at 3.5 cm is thus about 
23 mJy (Zurita et al. 2000) which is consistent with our previ- 
ous estimation of 6 mJy for the non-thermal radiocontinuum. 

Table |4] gives an estimate of the thermal and non-thermal 
contributions to the observed radio continuum emission, us- 
ing the Ha luminosities to estimate the thermal component. 
At 20 cm the radio continuum emission is mainly non-thermal 
(87 %), but at shorter wavelengths (6 and 3.5 cm) both contri- 
butions are comparable. 




Right Ascension (J2000) 



Fig. 8. Integrated '^CO(l-O) (Jy.beam^'km s"') emission in- 
tensity in the velocity range 20-100 km s^'overlaid on the 
Ha emission (greyscale). The contour levels are 0.5, 0.75, 
1, 2, 4, 6 and 10 Jy.beam^'km s"'(lcr rms noise ~ 0.02 
Jy.beam^'km s"'). 

4. Molecular gas 

4.1. Distribution 

The CO emission detected with the PdBI is strongly concen- 
trated toward the center of NGC 6764 (see Figure |2]i and 
reaches a maximum at a projected distance of r.'25 (~ 190 
pc) from the radio continuum peak. More than 98% of the to- 
tal molecular gas is concentrated within the inner 800 pc, and 
as much as 20% of the total '^CO( 1-0) emission is included 
within one beam. 

Apart from the central concentration, the CO emission is 
extended along the bar towards the western side ("W") and has 
a smaller extension to the eastern side. Assuming that the exter- 
nal spiral arms are trailing, this CO emission is on the leading 
side of the bar. 

The western CO feature, which extends to 10" from the 
center, is aligned with the primary bar The CO emission is 
associated with two Ha-emitting star forming regions (with 
5 X 10^^ and 3 x 10^' erg.s^^ Ha luminosities not corrected 
for extinction) which are among the brightest in NGC 6764, 
excluding the strong Ha emission at the tip of the bar (Rozas 
et al. 1996). Since NGC 6764 is highly inclined, the location of 
star formation and molecular gas in these regions is not com- 
pletely clear. 

Figure |7] shows the '^CO(l-O) channel maps at 10 km s ' 
velocity resolution relative to systemic. '^CO( 1-0) emission 
is detected in the individual channels at 2.5cr level from 
2210 to 2600 km s-'(LSR). Two i2cO(l-0) peaks are seen in 
the very center at negative and positive velocities. The first 
of these strong peaks is seen mainly at negative velocities 
(-110 to -20 km s '), and is separated from the primary CO 
peak by ~2". This secondary peak is not visible in the inte- 
grated '^CO(l-O) map. Extended weaker features are visible 
in the channel maps as well. One feature is seen north of the 
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19 08 17 16' 
R.A. 



Fig. 7. Channel maps for the '^CO(l-O) emission. The rms noise is 2 mJy beam The contours are at (3, 6, 12, 24, 48, 92) 
o". The synthesized beam is 2" x 1'.'6 in size. The rest velocity is 2400 km s"'. The cross is located at the phase center, i.e., at 
a/2000 =19:08:16.370 and 5y20oo=50:55:59.58. 



'^CO(l-O) peak at velocities -40 to 130 km s~' . Another more 
extended feature (see Figure [8]l is seen at velocities 20 to 100 
km s"'. The first feature is associated with the molecular out- 
flow, visible in the '^CO(l-O) integrated intensity image, and 
flaring up to a 8"(~ 1 kpc) height. Correcting for the disk in- 
clination, the molecular gas reaches 1.5 kpc above the plane of 
the galaxy. The velocities are consistent with motion perpen- 
dicular to the galaxy plane, taking into account the inclination 
and orientation of NGC 6764. 

The faint extended feature between 20 and 100 km s"' cov- 
ers an area of about 1 kpc^ between 5"and 10" radii (800 - 
1600 pc). At a velocity of 70 km s it reaches a total length 
of 12" (~1.9 kpc), and is oriented north to south. The velocity 
field shows few signs of outflow motion in this region, indicat- 
ing that this feature is in the plane of the galaxy. The western 



feature (W) has a large velocity extent from -80 km s~' up 
to 10 km s A CO spur is visible towards the north, and its 
location is spatially correlated with the boundary of the radio 
continuum superbubble, as seen at 20 cm. 

4.2. Kinematics 

4.2.1 . First and second velocity moments 

Figure |9] shows the beam-smeared map of the intensity- 
weighted '^CO(2-l) velocity field, derived from the first mo- 
ment of the '^CO(2-l) data cube, overlaid on the 20 cm ra- 
dio continuum emission. We determined the dynamical center 
of the gas motion from the '^CO(2-l) data, and compared it 
to the positions of the various emission peaks (see Table HJ. 
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Table 5. Position (J2000.0) of the peak emission in radio con- 
tinuum (20 cm), CO, and Ha (off-filter and the emission line). 
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HS06 find an offset of 0.55"between the VLA A-aiTay peak 
at 3.5 cm and the optical nucleus. They explain that difference 
by a higher extinction close to the nucleus, which is consistent 
with the molecular gas distribution. The dynamical center of 
CO emission is spatially coincident with the nucleus as seen 
in the B- and /-bands, and located about (X'S (~ 80 parsecs) 
from the radio continuum peak. Our observations are consis- 
tent with the global rotation of the galaxy disk, combined with 
streaming motions along the bar and/or a warping of the disk. 
The strongest deviation from the rotation pattern, the so-called 
"spider diagram," is seen on the eastern side along the molecu- 
lar outflow, which originates from the '^CO(l-O) peak, at 1"- 
2" to the NE of the radio continuum peak. The noise in the 
velocity field there is significant since the CO emission is cen- 
trally concentrated and the emission intensity is much lower 
beyond the peak. The position angle (P. A.) of the galaxy ma- 
jor axis, determined from the velocity field, 70° ±5, is slightly 
larger than the PA. from the orientation of the optical isophotes 
(61?5). The HI data of Wilcots et al. (2001) indicate a change in 
the PA. between the outer part (~ 58°) and the inner part of the 
disk (~ 76°). The difference is most likely due to the streaming 
motion near the center caused by the bar, as confirmed by the 
disturbed CO velocity field. 

Figure [TO] shows the beam-smeared map of the intensity- 
weighted '^CO(l-O) velocity dispersion in the center, derived 
from the second moment of the ''CO(l-O) data cube, overlaid 
with the '^CO(l-O) velocity field. The central velocity disper- 
sion (Tv - 200 km s is much higher than the internal veloc- 
ity dispersion of a molecular cloud (typically a few km s"'; see 
e.g. Scoville et al. 1987). More likely, it results from the beam- 
smeared velocity gradient caused by a steeply rising rotation 
curve, together with a contribution from streaming motions. 

Within the astrometric accuracy the peak of the velocity 
dispersion is located at the position of the dynamical center. At 
the position of the molecular outflow the velocity dispersion 
has values up to 150 kms"', whereas along the bar, at radii 
larger than 5" from the nucleus, the molecular gas velocity dis- 
persion does not exceed 50 km s"' . Curiously, the velocity dis- 
persion is asymmetric in the east-west direction, with larger 
values on the eastern side. 




19''08"'17.0= 



165= 16-0= 
Right Ascension (J2000) 



Fig. 9. Isovelocity contours of the '^CO(2-l) emission over- 
laid on a color image of the 20 cm radio continuum emission. 
Black contour lines indicate negative velocities, and are sep- 
arated by steps of 10 km s White contour lines imply pos- 
itive velocities with respect to the reference velocity, and are 
also separated by steps of 10 km s"'. The reference velocity is 
2400 km s-'. 



+ 50^56'10" — 




191^08^17. 



Fig. 10. Velocity dispersion of the '^CO( 1-0) emission in 
greyscale, overlaid by the '^CO(l-O) velocity field. The white 
and black lines are for positive and negative velocity relative to 
the reference velocity at 2400 km s ' . 



4.2.2. Lindblad resonances 

The rotation curve has been computed from our 
'2CO(1-O)data and the HI data from Wilcots et al. (2001), 
using the task INSPECTOR in Gipsy. Although the HI spatial 
resolution is poorer than in CO, it provides data points at 
large radii. The rotation curve becomes flat at a value of 250 
km s ' at ~ 12 kpc from the center Under the assumption 
of the epicyclic approximation, the azimuthal frequency (Q.) 
and the radial frequency (k) have been computed from the 
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Table 6. Dynamical resonances as estimated from a rough 
epicyclic model for NGC 6764. Q.i,„r-22 km s ' kpc '. See the 
text for a detailed description. 



Resonance 


n 


Radius 




(Q.l,ar = n + nX k) 




ilLR 


-1/2 


6±6" 


oILR 


-1/2 


10±10" 


iUHR 


-1/4 


24±5" 


CR 





40±12" 


oUHR 


+ 1/4 


100±10 



rotation curve. The bar pattern speed is estimated from the 
rotation curve and the length of the bar semi-major axis (~6.5 
kpc). Several studies (e.g., Athanassoula 1992, Combes & 
Elmegreen 1993, Sempere et al. 1995) have shown that the 
ratio of the bar corotation radius to the bar semi-major axis 
length is in the range 1-1.5. The determination of the bar 
pattern speed is not strongly dependent on the exact value of 
this ratio since the frequency Q is flat in that domain. The bar 
pattern speed Q^ar is ~22 km s"' kpc"' with a corotation (CR) 
at a radius of 40" (6.5 kpc), using a value of 1 for the ratio of 
the radius of the CR to the length of the bar semi-major axis. 
The inner Lindblad resonances (ILRs), if they exist, are the 
loci where the bar pattern speed equals Q. - k/2 (Figure fTTT i. 
Thus, the epicyclic approximation suggests that two ILRs 
are present near the center. The rough radii for the inner and 
outer ILRs (ilLR and oILR) are 6" and 10" (1 and 1.7 kpc). 
4:1 resonances, or (inner/outer) ultraharmonic resonances 
(i/oUHR), with Q.bar ^ Q. + k/4, are found at radii of 24" (3.9 
kpc; iUHR) and 100" (16.5 kpc; oUHR). At these resonances 
the star orbits are periodic and form rings (e.g, Byrd et al. 
1998). Contrary to the UHR location, the locations of the ILRs 
strongly depend on the bar pattern speed. Table |6] summarizes 
the locations of the resonances found with the application of 
the epicyclic approximation, valid for asymptotically weak 
bars. In the presence of an even moderately strong bar, an 
accurate determination of the resonance positions requires a 
full dynamical model, together with orbit analysis, which is 
beyond the scope of this paper. 

4.2.3. Non-circular motions 

Figure [12] shows the position-velocity (PV) diagrams of the 
'^CO(l-0)and '^CO(2-l) emission along the major and mi- 
nor axes (RA. = 62° and RA. = 152°) of NGC 6764, along 
with a PV diagram at PA. = 90° for comparison with the 
'^CO(l-O) observations made with the Nobeyama Millimeter 
Array (Eckart et al. 1996). As already shown by Eckart et 
al. (1996), the PV diagrams imply an asymmetric distribu- 
tion of molecular gas near the nucleus. The CO intensity is 
higher for gas at redshifted velocities, with respect to the sys- 
temic velocity, i.e., on the eastern side. The PV diagram at 
P.A. = 90° in the '^CO( 1-0) emission shows strong stream- 
ing motion at x=-4", which is located slightly inside the 
ilLR (r ~ 6"). The '2cO(2-l)PV diagram along the minor 
axis reveals large streaming motions at x = l"-2", with a 
velocity near -200 kms Similarly, the '2CO(1-0)PV di- 




R (kpc) 

Fig. 11. Azimuthal frequency Q (stars) estimated from the ro- 
tation curve (HI-hCO data). The first three data points are from 
our new CO (1-0) observations. The data points at larger radii 
are from HI observations (Wilcots et al. 2001). The UHR (dia- 
mond) and ILR (cross) locations have been estimated by set- 
ting the bar corotation radius to the semi-major axis radius 
of the bar (~40"). The bar pattern speed is estimated to be 
22 km s"' kpc"^ 



agram along the minor axis shows an area of lower density 
gas on the eastern side at radii < 5", with streaming motions 
at negative velocities. The spur visible between x-positions 
1" and 6" in the velocity range [-30,80]km s"' (relative to 
2400 km s"') is located in the molecular outflow, and is con- 
sistent with the '^CO(l-O) spectrum integrated over the out- 
flow area. Assuming a trailing spiral structure, the northern 
part of NGC 6764 is closer to us, which implies largely red- 
shifted (positive) velocities for the radial outflow in the north- 
ern side, as is seen in the PV diagram and in the integrated 
'2CO(2-l) spectrum. 

4.3. CO line ratio 

The line ratio for the integrated CO intensity be- 
tween the '2CO(2-l) and ^^CO(l-O) transitions, 

= /(>2cQ(2_i))//(i2co(l-0)), was computed matching 
the u-v coverage of both transitions and using the same spatial 
resolution and taking only the '^CO(l-O) and '^CO(2-l) data 
above a threshold of 3cr. The result is shown in Figure [T3] 
The maximum % is 0.98, in agreement with that estimated by 
Eckart et al. (1991) from single dish observations. Eckart et 
al. (1991) interpreted the global ratio to arise from moderately 
warm, dense, and optically thick gas. 

To estimate the uncertainty in the line ratio introduced by 
the missing flux, we compared the CO line ratio map derived 
from the single-dish data with the ratio derived from our new 
PdBI observations. We derive larger values of the ratio at large 
radii, owing to the lower spatial resolution (~ 20") of Eckart 
et al. (1991). Their maximum value at the center is lower than 
ours by ~ 10%. The maximum line ratio !R,„o.v in each channel 
(10 km s ' width) is shown in Figure [14] A bimodal feature is 
clearly seen, with peaks towards the blue and red sides. The 
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Fig. 12. Position-velocity diagram of the ^^CO(l-O) (left) and ^^€0(2-1) (right) emission in the center of NGC 6764. The zero- 
offset is at the CO peak. The major axis P. A. is 62°, and the PA. of 90° has been included to compare with Fig. 1 of Eckart et al. 
(1996). The units are in Jy beam 



blue side with a peak at -140 km s" is associated with molec- 
ular gas closer to the radio continuum peak, as shown in the 
CO line ratio channel maps in Figure [15] At the peak the ratio 
reaches values of 2-2.5. The red side of the peak at 140 km s"' 
has a maximum value of about 1.3. This high CO line ratio in 
fact applies to the bulk of the molecular gas, as is shown by the 
CO line ratio channel maps below. 

The '2CO(2-1)/'2CO(1-O)line ratio has been computed 
for each channel (10 km s 'width). Figure [TSl shows the CO 



line ratio from a channel at -180 km s ' up to a channel at 200 
km s"' (the systemic velocity of NGC 6764 is 2400 km s"', 
corresponding to a zero velocity in the channel maps). From 
-180 to -40 km s ' the peak of the CO line ratio is spatially 
coincident with the nucleus (or the radio continuum peak, as in- 
dicated by a cross in the figure). The peak at -140 km s"' seen 
in Figure[T4]is likely to be an edge effect. The "true" maximum 
is reached at -130 km s"', where the ratio reaches a value of 
1 .4 at the exact position of the radio continuum peak (to within 
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Fig. 15. '^CO(2-l) to '^CO(l-O) line ratio for different velocities, as indicated in each panel, with a channel width of 10 km s ' . 
The first contour is at 0.4 and the contours increase in steps of 0.1. The radio continuum peak is shown by a cross. 



a pixel width). The western CO emission feature has a mean 
CO line ratio of 0.3, which is likely to be an underestimate be- 
cause of the filtering out of the extended emission. 

The high CO line ratios for the blueshifted gas in the nu- 
cleus may be associated with outflowing molecular gas pow- 



ered by the AGN or the starburst. HS06 reported a possible 
blueshifted HI absorption at ~ 120km s"' against the compact 
core which may be related to the molecular gas with such high 
CO line ratios. We note that the radio continuum map (Fig.[T) 
has hints of a nuclear radio jet towards the south, likely associ- 
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Right Ascension (J2000) 



Fig. 13. Integrated intensity CO(l-O) line (greyscale) over- 
laid by contours showing the '^CO(2-l) to '^CO(l-O) line 
ratio, integrated over all channels. The contours are 0.4, 0.5, 
0.6, 0.7, 0.8,0.9 and 0.95. 
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Fig. 14. Maximum i2cO(2-l) to '^COCI-O) line ratio ^R„„., 
computed in the center of NGC 6764 for each '^CO(l-O) chan- 
nel (10 km s' width). 

ated with blue-shifted motions there as expected if the orienta- 
tion of NGC 6764 with respect to the line of sight is taken into 
account. A redshifted counterpart with such a high line ratio is 
not detected. 

4.4. mass 

Adopting a standard conversion factor of Xqq = 2.3 x 
10^°cm"2 (Kkms"')"' (e.g.. Strong etal. 1988), the total mass 
of molecular hydrogen (H2) is: 

Mh, = 9 X 10^ S CO (2) 

where S co is the integrated CO flux density in Jy km s ' and D 
is the distance in Mpc. The integrated '^CO(l-O) flux density 
from our new map gives a total molecular hydrogen mass of 
1.0 X 10^ Mq. It does not take into account the 50 % of flux 
filtered out by the interferometer observations. 



In the Milky Way's disk the Xqq conversion factor derived 
from observations has been found to match the theoretical con- 
version factor quite well (e.g., Radford et al. 1991). However, 
there is growing evidence that the conversion factor is overes- 
timated by a factor of 3-10 in the central part of galaxies, es- 
pecially in those with nuclear starbursts (Downes & Solomon 
1998). The physical reasons for this discrepancy have been 
discussed in several papers, and include the metallicity gradi- 
ent (e.g., Dahmen et al. 1998), broadening of the CO line by 
the bulge potential (e.g., Stacy et al. 1989, Mauersberger et al. 
1996), and changes in the physical conditions of the circumnu- 
clear gas (e.g., its optical depth, temperature, and density). In 
NGC 253, a nearby starburst galaxy, Mauersberger et al. (1996) 
found a conversion factor ~ 7 times lower than X^'^ in the cen- 
tral part of the galaxy. They argued that the central potential 
well is broadening the CO lines, and its emission could be par- 
tially optically thin. The conversion factor has been found to 
be 3 to 10 times lower in the galactic bulges compared to their 
disks (Sodroski et al. 1995). Recently, WeiB et al. (2001) ana- 
lyzed the variations of X(^q in the starburst of M82 in detail. 
They found that " X(^q is a function of the intrinsic gas proper- 
ties, which strongly depend on environmental effects" (massive 
star formation). This leads to a variation of X by a factor of 5- 
10, compared to the Galactic value of X(^q = 1.6 x 10^" cm"^ 
Kkms-i. 

The dust mass computed from the IRAS fluxes at 60 and 
100 //m gives a total mass of M.d-\A x 10^ Al© with a mean 
temperature of 36.5 K. It is known from FIR observations with 
ISO that the IRAS data are missing a cold dust component that 
emits beyond 150 ^va and has a temperature of Id < 20 K, 
(Alton et al. 1998, Trewhella et al. 2000). This component can 
be 2-10 times more massive than the warm dust component. 

In NGC 6764 the gas-to-dust mass ratio, using the standard 
CO-to-H2 conversion factor and correcting for the missing flux 
in the interferometer map, leads to Mjj^/M,;,,,, ~ 1500 + 225. 
This factor is much larger than the typical value found in nor- 
mal spiral galaxies (an average ratio of ~ 700; Leon et al. 1998) 
using IRAS-derived dust masses. This indicates that the molec- 
ular gas mass obtained from the '^CO(l-O) luminosity and the 
standard conversion factor has most likely been overestimated. 

The high gas-to-dust mass ratio means that the H2 mass has 
probably been overestimated by a factor of 2-3 in the center 
of NGC 6764, when using the standard conversion factor X~ 
X^Q. This overestimate may be even more severe in the nucleus 
because of the intense starburst. Thus we estimate the total H2 
mass to be 7xlO^Al0, taking into account the filtered out flux 
in the '^CO(l-O) interferometric map. 

5. Outflow 

The basic morphological parameters for the outflow from the 
20 cm radio continuum observations are: a projected base of ra- 
dius 3.5"(540 pc) and a projected height of 7"(1080 pc) above 
the galaxy plane. The mass-load in the outflow is estimated 
from the spatial distribution of the '^CO(l-O) emission. The 
'^CO(l-O) emission was integrated in the area shown in Fig.|2] 
only for the northern part of the outflow. The area was chosen in 
order to avoid the contamination from the molecular gas emis- 
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sion in the disk, given the incUnation of the galaxy and the spa- 
tial resolution of the observations. The '^CO(l-O) integrated 
intensity in the out-of-plane outflow is 1.4 Jy km s"'. 

The observed energy of the different gas phases in the out- 
flow is an indicator of the physical processes responsible for 
the outflow. Assuming that the bulk of the molecular gas is 
perpendicular to the disk plane, the expansion velocity Nexp 
is estimated from the spur in '^CO(l-O) seen between -30 and 
80 km s"^on the PV diagram (see Fig. ). We take a projected 
expansion velocity of 25 km s"\ which gives an expansion ve- 
locity of 53 km s"' . This value is quite uncertain since we know 
already that some molecular gas is reaching an expansion ve- 
locity of 170 km s"'. Applying a correction of 1/3 to the CO- 
to-H2 conversion factor (see §4.4) the total molecular gas mass 
in the outflow is 4.3 x 10* Mq. Thus the molecular gas kinetic 
energy Efa„ in the outflow is = \Mgas * V^^p = 2.4 x 10^^ 
ergs. 

For a starburst with a SFR=1 Alg/yr, different star forma- 
tion models have been computed (instantaneous or continu- 
ous). The total energy released by SNs and stellar wind at the 
age of the recent NGC 6764 starburst (~ 5 Myr) is in the range 
of 10^^ - 10^^ ergs. Thus there is no need to invoke AGN or 
other mechanisms, such as hypernovae or tidal interaction, to 
drive the outflow, even if the AGN seems to have an effect 
on the molecular gas physical conditions (see Sect. 4.3) and 
.the data cannot exclude the AGN from being the major energy 
source either. 

To further study the energy budget, we computed the total 
thermal energy in the hot X-ray gas, which was determined 
by the HRI ROSAT instrument. The X-ray emission is com- 
posed by two components: one varying on the timescale of sev- 
eral days and a more diffuse one which extends up to 15" from 
the nucleus, roughly along the minor axis (Schinnerer et al. 
2000). We make use of the X-ray parameters derived by Eckart 
et al. (1996) from the /?05Ar observations in the 0.1-2.4 keV 
band. They estimated an electron density He ~ 0.2/7''^cm"^, 
where fs is the volume filling factor, the volume of the extended 
X-ray gas as V ~ 7 x 10''"' cm-', and the temperature of the X- 
ray gas as T ~ 4 x 10* K. The total thermal energy is then 
given by ~ in^kTV/s, i.e. ~ 2.3 x 10^ fj'^ ergs. We 
use /, = 0.01, as was found for the starburst galaxy NGC 3256 
(Heckman et al. 1990; Moran et al. 1999). Thus Ejc is estimated 
to be ~ 2.3 x 10^^ ergs. The current starburst in NGC 6764 ap- 
pears to be capable of releasing sufficient energy to power the 
hot X-ray gas bubble. 

To estimate the distribution of the radio continuum emis- 
sion at 20 cm in the outflow bubble, the nuclear component 
has been subtracted by fitting a Gaussian to the nuclear flux 
above a threshold of 4 mJy beam"' . The total flux remaining in 
the center-subtracted map is 63 mJy, or about 60% of the total 
seen in the VLA observation. There is an offset of the centroid 
of the outflow emission relative to the position of the nuclear 
radio continuum peak, and relative to the disk plane. This off- 
set can be due to various reasons: a) the starburst is located in 
the western part of the nucleus, and sUghtly below the plane, 
as suggested by the lower gas column density south of the nu- 
cleus, b) there is a possible contribution from a radio jet in the 
southern part in addition to the starburst-induced outflow. 



6. Discussion 

We compare the CO and Ha observations to the different hy- 
drodynamical models of Tomisaka & Ikeuchi (1988) performed 
for a smooth gas component with various scenarios of starburst. 
The case D of Tomisaka & Ikeuchi (1988) has a starburst with 
a supernova rate comparable to the one found in NGC 6764 
(lO^^yr"': the bubble reaches z~1.2 kpc after 9 Myr, which 
is difficult to conciliate with the W-R star age and population. 
In the case C, with an ambient density of n « 100 cm"^, the 
height of the bubble reaches 1 kpc after 7 Myr Moreover they 
performed a simulation with a short "explosion" which drives 
the bubble height to z « 1.3 kpc as fast as in 1 Myr. From all 
these different models, we can conclude that the recent starburst 
has been intense and short 3-7 Myr ago, in complete agreement 
with the stellar population analysis of Schinnerer et al. (2000). 
From the comparison with these simulations we can conclude 
that the hot gas in NGC 6764 is still in a confined bubble, and 
has not yet escaped to the halo. 

The NSB in NGC 6764 is very similar to the NSB in 
NGC 2782 (Jogee et al. 1998, 1999). However, NGC 2782 has 
no nuclear activity, apart from the starburst. The CO spurs in 
NGC 2782 have a total molecular gas mass of 2 X 10^ Al©, a 
total molecular gas kinetic energy of ~ 8 x 10^^ ergs, and a 
thermal energy of ~ 2 x 10^*/*^. These numbers are larger by 
factors of 5-10 compared to the energy released in NGC 6764. 
Nevertheless the size of the outflows is comparable. This is 
probably due to the lower ISM density in NGC 6764 than in 
NGC 2782. 

Observations (Seaquist & Clark 2001, Walter et al. 2002) 
of '^CO(3-2) and '^CO(l-O) emission in the prototypical star- 
burst galaxy M82 indicate a large amount of molecular gas 
(> 3 x 10*^ Al© in the outflow with a maximum outflow velocity 
of 230 km s"'and an estimated kinetic energy of ~ IQ^^ergs. 
Similarly to our estimate, WeiB et al. (2001) found a CO-to-H2 
conversion 3 times lower than the Galactic conversion factor. 
The height of the molecular gas outflow is about 1 kpc above 
the plane with an estimated starburst age of 10 Myr. The com- 
parison with NGC 6764 indicates that in this galaxy the build- 
up of the molecular outflow was faster by a factor 2-3. The 
main difference is the molecular gas mass involved in the out- 
flow, which is an order of magnitude larger in M82 than in 
NGC 6764. 

HS06 present a compilation of 10 galaxies with non- 
thermal bubbles that shows that all are associated with an AGN. 
They link the origin of these bubbles to the presence of an 
AGN. We note that 6 galaxies of their sample (60 %) have a 
FIR luminosity larger than the one of NGC 6764, if we adopt 
their luminosities. The presence of an AGN may not be suffi- 
cient to create these bubbles. It may indicate that the associa- 
tion of a starburst, present or recent, with an AGN could be the 
key to generate non-thermal radio bubbles. The creation of such 
feature would depend on a short time scale combination of the 
two processes: the starburst and the AGN. Nevertheless we note 
that the presence of radio AGN is very low in spiral galaxies (J. 
Sabater, private communication) which makes the correlation 
between radio- AGN and non-thermal bubble very significant. 
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Thus a statistical study with a large sample of galaxy with non- 
thermal bubbles would help to disentangle their origins. 

NGC 6764, hke other composite galaxies with a nuclear 
starburst and an AGN, shows that the molecular gas evolution 
in the very center is mainly driven by the starburst, especially in 
the molecular outflow. Even if the AGN has an active role in the 
central evolution, as shown by the gas chemistry in the Seyfert 
2 galaxy NGC 1068 (Usero et al. 2004), the only possible in- 
fluences of the AGN on the molecular gas, at the resolution 
of the observations of the CO emission in NGC 6764, appears 
to be the high CO line ratio in the very center, and likely the 
high HCN/CO ratio already found by Kohno et al. (2001). No 
clear dynamic signature is observed. The AGN are able to in- 
ject vast amounts of energy into their host galaxies and should 
have a disruptive influence on the molecular gas in the central 
regions. 

7. Summary 

We have found the following results from our high-resolution 
mapping of the '^COCl-O) and ^^00(2-1) milhmetre lines 
and 3.5 and 20 cm radio continuum emission in the circum- 
nuclear area of the NSB galaxy NGC 6764. 

1. The bulk of the '^CO emission is concentrated in the 
very center of NGC 6764 and is offset from the 20 cm 
radio continuum peak by ~ I'.' 2 (with respect to the 
'^CO(l-O) emission peak). A more extended component is 
detected along the bar on the western side of the nucleus. 

2. Non-thermal radio continuum emission at 3.5 and 20 cm is 
observed above and below the galactic plane in a nuclear 
bubble-like outflow, as already shown by HS06 and Laine 
et al. (2006). This outflow is detected out to z=1.3 kpc 
above the galactic plane (deprojected). ^^CO(l-O) emission 
is detected at the NE boundary of the outflow. About 
4.3 X 10* Mq of molecular gas is driven by the outflow 
out of the plane of the galaxy. 

3. The total molecular gas mass is 6.7 x 10^ Al© using a CO- 
to-H2 conversion factor equal to 1/3 of the Galactic value. 
The use of such a value is motivated by the large gas-to- 
dust mass ratio and previous studies of molecular gas in 
starburst galaxies. 

4. The i2CO(2-l) / line ratio was computed in 
each channel (10 km s^'velocity width). It reaches a maxi- 
mum of 1-2 at the nucleus, and is likely associated with the 
nuclear AGN (it coincides with the 20 cm radio continuum 
peak). In the bulk of the molecular gas the maximum ra- 
tio is about 1-1.3, consistent with warm (T4,„ > 20 K) and 
dense gas (njj^ > 2 X lO'^cm"^), as found by Eckart et al. 
(1991). 

5. An important fraction of the radio continuum at 3.5 cm 
is thermal. The outflow morphology and the ISM proper- 
ties indicate an age of ~ 3-7 Myr for the recent starburst, 
by comparing the observations to hydrodynamical simula- 
tions. It appears that the outflow is still confined and not 
freely expanding. A "hot-spot" of thermal gas with a flat 
spectral index ( 0.2) is located near the northern bound- 
ary of the bubble/outflow. 



6. The kinetic and thermal energies of the molecular and hot 
gas components in the bubble can be fully accounted for by 
the energy released from the NSB, even if the data cannot 
exclude the AGN from being the major energy source. 

Acknowledgements. We thank an anonymous referee for his care- 
ful reading and very detailed report which helped to improve this 
paper significantly. This research has made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. During part of 
this research S. Leon was supported by a Marie Curie Individual 
Fellowship contract HPMF-CT-20002-01734 from the European 
Union. We would like to thank David Prayer and Eric Wilcots for pro- 
viding their Ha and HI data, respectively. 

References 

Alton, P B., Bianchi, S., Rand, R. J., Xilouris, E. M., Davies, J. I., 

Trewhella, M. 1998, ApJ, 507, 125 
Athanassoula, E. 1992, MNRAS, 259, 345 

Baum, S. A., O'Dea, C. P., Dallacassa, D., de Bruyn, A. G., Pedlar, A. 
1993 ApJ, 419, 553 

Byrd, G.G., Ousley, D., dalla Piazza, C. 1998, MNRAS, 298, 78 

Caplan, H., Deharveng, L. 1986, A&A. 155, 297 

Colina, L., Arribas, S., Borne, K.D. 1999, ApJ, 527, L13 

Combes, F. 2001, 'The Central Kiloparsec of Starbursts and AGN: 
The La Palma Connection', ASP Vol. 249. Eds. Knapen, J.H., 
Beckman, J.E., Shlosman,I, Mahoney, T.J, p. 475 

Combes, E, Gerin, M. 1985, A&A, 150, 327 

Combes, E, Elmegreen, B.G. 1993, A&A, 271, 391 

Condon, J. J., Yin, Q. E 1990, ApJ, 357, 97 

Dahmen, G., Huttemeister, S., Wilson, T. L., Mauersberger, R. 1998, 

A&A, 31, 959 
Dopita, M.A. 1985, ApJ, 295, 5 
Downes, D., Solomon, R M. 1998, ApJ, 507, 615 
Eckart, A., Cameron, M., Boiler, Th., Krabbe, A., Bhetz, M., Nakai, 

N., Wagner, S. J., Sternberg, A. 1996, ApJ, 472, 588 
Eckart, A., Cameron, M., Jackson, J. M., Genzel, R., Harris, A. I., 

Wild, W, Zinnecker, H. 1991, ApJ, 372, 67 
Franceschini, A., Silva, L., Fasano, G., Granato, L., Bressan, A. 

Arnouts, S., Danese, L. 1998, ApJ, 506, 600 
Hota, A, Saikia, D.J. 2006, MNRAS, 371, 945 (HS06) 
Heckman, T. M., Armus, L., Miley, G. K. 1990, ApJS, 74, 833 
Hummel, E., van Gorkom, J. H., Kotanyi, C. G. 1983, ApJ, 267, 5 
Hummel, E., Beck, R, Dettmar, R.-J. 1991, A&AS, 87, 309 
Hill, T. L., Heisler, C. A., Norris, R. P, Reynolds, J. E.; Hunstead, R. 

W. 2001, AJ, 121, 128 
Irwin, J. A., Sofue, Y. 1996, ApJ, 464, 738 
Jogee, S. Kenney, J. D. R, Smith, B. J. 1998, ApJ, 494, 185 
Jogee, S., Kenney, J. D. P, Smith, B. J. 1999, ApL 526, 665 
Kenney, J. D. P, Carlstrom, J. E., Young, J. S. 1993, ApJ, 418, 687 
Kennicutt, R. C, Jr. 1998, ApJ, 498, 541 

Kohno, K. et al. 2001, "The Central Kiloparsec of Starbursts and 
AGN: The La Palma Connection", ASP Conference Proceedings 
Vol. 249. Eds. Knapen, J. H., J. E. Beckman, J.E., Shlosman, I. and 
Mahoney, T.J. San Francisco: Astronomical Society of the Pacific, 
2001, p. 672. 

Kunth, D. Legrand, F. Tenorio-Tagle, G. SiUch, S. Mas-Hesse, J. M. 

Cervifio, M. 2002, Ap&SS, 281, 261 
Krolik, J.H., Begelman, M.C. 1986, 308, L55 

Laine, S., Kotilainen, J. K., Reunanen, J., Ryder, S. D., Beck, R. 2006, 
AJ, 131, 701 



Leon et al.: Molecular gas in NGC 6764 



15 



Laine, S., Shlosman, I., Knapen, J. H., Peletier, R. F. 2002, ApJ, 567, 

97 

Launhardt, R., Zylka, R., Mezger, R G. 2002, A&A, 384, 112 

Lee, S.-W., Irwin, J. A., Dettmar, R.-J., Cunningham, C. T., GoUa, G., 

Wang, Q. D. 2001, A&A, 377, 759 
Leon, S., Combes, R, Menon, T. K. 1998, A&A, 330, 37 
Leon, S., Meylan, G., Combes, F. 2000. A&A, 359, 907 
Mauersberger, R., Henkel, C, Wielebinski, R., Wiklind, T., Reuter, 

H.-R 1996, A&A, 305,.421 
Mihos, J. C, Hemquist, L. 1994, ApJ, 425, 13 
Moran, E. C, Lehnert, M. D., Helfand, D. J. 1999, ApJ, 526, 649 
Norman, C. A., Ikeuchi, S. 1989, ApJ, 345, 372 
Osterbrock, D. E., Cohen, R. E. 1982, ApJ, 261, 640 
Radford, S. J. E., Downes, D., Solomon, R M. 1991, ApJ, 368L, 15 
Rozas, M.; Beckman, J. E.; Knapen, J. H. 1996, A&A, 307, 735 
Sakamoto, K. Ho, R T. P., lono, D., Keto, E. R.; Mao, R.-Q., 

Matsushita, S., Peck, A. B., Wiedner, M. C; Wilner, D. J., Zhao, 

J.-H. 2006, ApJ, 636, 685 
Sanders, D. B., Mirabel, I. F. 1985, ApJ, 298, 31 
Scalo, J., Chappell, D. 1999, MNRAS, 310, 1 

Scoville, N. Z., Yun, M. S., Sanders, D. B., Clemens, D. R, Waller, W. 

H. 1987, ApJS, 63, 821 
Schinnerer, E., Eckart, A., Boiler, Th. 2000, ApJ, 545, 205 
Seaquist, E. R., Clark, J. 2001, ApJ, 552, 133 

Sempere, M.J., Garcia-Burillo, S., Combes, F., Knapen, J. H. 1995, 

A&A, 296, 45 
Sodroski T.J. et al. 1995, ApJ, 452, 262 

Stacy, J. G., Bitran, M. E., Dame, T. M., Thaddeus, R 1989 The Center 
of the Galaxy: Proceedings of the 136th Symposium of lAU, Los 
Angeles, U.S.A., Ed M. Morris, Kluwer Academic PubUshers, 
Dordrecht, p. 157 

Strickland, D. K., Ponman, T. J., Stevens, I. R. 1997, A&A, 320, 378 
Strong, A. W., Bloemen, J. B. G. M., Dame, T. M., Grenier, L A., 

Hermsen, W., 

Taylor, D., Tyson, Axon, D.J. 1992, MNRAS, 255, 351 

Tomisaka, K., Ikeuchi, S. 1988, ApJ, 330, 695 

Trewhella, M., Davies, J. L, Alton, P. B., Bianchi, S., Madore, B. F. 

2000, ApJ, 543, 153 
Usero, A., Garca-Burillo, S., Fuente, A., Martn-Pintado, J.; Rodrguez- 

Femndez, N. J. 2004 A&A, 419, 897 
Walter, F.,Weiss, A., Scoville, N. Z. 2002, ApJ, 580, 21 
WeiB, A., Neininger, N., Huttemeister, S., Klein, U 2001, A&A, 365, 

571 

White, R. L., Becker, R. H. 1992, ApJS, 79, 331 

Zurita, A., Rozas, M., Beckman, J. E. 2000, A&A, 363, 9 



